Fortschritte der Physik, 2 March 2013 



Preprint: LA-UR- 12-20346 

Non-Hermitian approach for modeling of noise-assisted quan- 
tum electron transfer in photosynthetic complexes 

Alexander L Nesterov 1 *, Gennady P. Berman 2 **, and Alan R. Bishop 3 *** 

1 Departamento de Ffsica, CUCEI, Universidad de Guadalajara, Av. Revolution 1500, Guadalajara, CP 
44420, Jalisco, Mexico 

2 Theoretical Division, MS-B213, Los Alamos National Laboratory, Los Alamos, NM 87544, USA 
^Sj ", 3 STE, MS-A127, Los Alamos National Laboratory, Los Alamos, NM, 87544, USA 

7— i ; 

O ! Received XXXX, revised XXXX, accepted XXXX 

(N ■ Published online XXXX 

a: 



< 



43 

Or 



o 



Key words non-Hermitian Hamiltonian, reaction center, electron transfer, noise 



We model the quantum electron transfer (ET) in the photosynthetic reaction center (RC), using a non- 
Hermitian Hamiltonian approach. Our model includes (i) two protein cofactors, donor and acceptor, with 
discrete energy levels and (ii) a third protein pigment (sink) which has a continuous energy spectrum. Inter- 
actions are introduced between the donor and acceptor, and between the acceptor and the sink, with noise 
acting between the donor and acceptor. The noise is considered classically (as an external random force), 
and it is described by an ensemble of two-level systems (random fluctuators). Each fluctuator has two in- 
dependent parameters, an amplitude and a switching rate. We represent the noise by a set of fluctuators 
with fitting parameters (boundaries of switching rates), which allows us to build a desired spectral density 
of noise in a wide range of frequencies. We analyze the quantum dynamics and the efficiency of the ET as 
a function of (i) the energy gap between the donor and acceptor, (ii) the strength of the interaction with the 
continuum, and (iii) noise parameters. As an example, numerical results are presented for the ET through 
the active pathway in a quinone-type photosystem II RC. 
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Nature has evolved photosynthetic organisms to be extremely complex bio-engines that capture visible 
, light in their peripheral light-harvesting complexes (LHCs) and transfer excited-state energy (as excitons) 

through the proximal LHC of photosystem II (PSII) and photosystem I (PSI) to the RCs. The primary 
charge separation occurs in the RC (which works as a battery), leading to the formation of an electro- 
chemical gradient 0J2). During the past two decades, crystallographic structures for many photosynthetic 
, complexes (PCs), including the LHCs and RCs, have been determined to a resolution of 2.5-3 A JSpO. 

(See also references therein.) 

Like all engines, PCs operate in a thermal environment at ambient temperature and in the presence of 
external "classical" sources of noise ll8tiT4l. In spite of this, recent experiments based on two-dimensional 
laser-pulse femtosecond photon echo spectroscopy revealed a long-lived exciton-electron quantum coher- 
ence in PCs such as the Fenna-Matthews-Oslov (FMO) and marine algae lfT5l - fl7l . Mainly, this occurs 
because the dynamics of the ET is so rapid (some picoseconds) that the thermal fluctuations and external 
noise are unable to significantly destroy quantum coherence. Consequently, the exciton/electron dynam- 
ics in LHCs-RCs must be described using quantum-mechanical methods 11 1 844221 . (See also references 
therein.) An important consequence of this is the high ET efficiency of the peripheral antennae complexes 
(close to 100 %). 
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As an example, Fleming and colleagues [18] have modeled quantum coherence effects in the bacterial 
FMO LHC by (i) using a tight-binding model (TBM) for exciton dynamics and (ii) introducing an empirical 
thermal relaxation function having an exponential form, in order to describe the high efficiency of exciton 
energy transport. Usually, in the TBM the exciton/electron ET dynamics in LHCs-RCs is described in the 
single exciton/electron approximation (due to limited sunlight intensities), with N ( N = 7 for the FMO in 
lfT8l ) being the total number of discrete pigments/sites. (Note that more complicated models which account 
for exciton and charged states can also be used El .) In this case, each pigment, n (n = 1, . . . , N), is 
represented by a two-level system with states |0 rl ) (unoccupied) and |1„) (occupied). The total Hamiltonian 
is H to t — H e + H p h + H e i- p h |fT8l . The first term is the Hamiltonian of exciton/electron states of the 
pigments in the site representation: H e = J2i E n \n)(n\ + ~Y^ n ^ n V mn \m) {n\, where E n is the site 
energy, and V mn denotes the coupling between the n-th and m-th pigments. The term H p h describes the 
thermal phonons provided by the protein environment, and the third term describes the interaction between 
pigments and the thermal phonons. It was numerically demonstrated in |18|, that in the FMO complex, 
quantum coherent ET is an adequate way to describe the energy transport dynamics. 

Usually, there are two different approaches which are used to describe the influence of the protein 
environment on the ET. One is based on the thermal environment J9)- In this case, the environment 
acts self-consistently on the electron system and, in combination with the transition amplitudes between 
sites/pigments, provides the ET rates between the sites and the Gibbs equilibrium state for the LHC-RC 
subsystem. The other approach is based on considering an external "classical" noise ED provided by the 
protein vibrations. This approach results in a transfer rate for the electron, but does not lead to Gibbs equi- 
librium states. The choice of approach depends on the concrete experimental situation which the theoretical 
model is intended to describe. 

In this paper we use the second approach, modeling the noise by an ensemble of fluctuators lfl4l . To 
simplify our description, we introduce a set of fluctuators with fitting parameters (boundaries of switching 
rates between relatively slow and fast fluctuators), which allows us to build a desired spectral density 
of noise in a wide range of frequencies. In particular, the spectral density of noise, used in this paper, 
includes the components of white noise, 1/ / noise, and high-frequency noise. We demonstrated in [ 14] that 
this approach successfully described the experiments El on the quantum dynamics of superconducting 
qubits. Here we consider the simplest model of ET in a quinone-type active pathway of the PSII RC. Our 
model includes two protein cofactors (donor and acceptor) with discrete energy levels, with the acceptor 
being embedded in a third protein pigment (sink) that has a continuous energy spectrum. In EH an 
additional sink reservoir was empirically introduced in order to describe the high ET efficiency in the FMO 
complex. A sink reservoir was also introduced phenomenologically in (2T| to describe the asymmetry of 
two branches of the ET in the photosynthetic PSII RC, and in j24j to describe the dynamics of excitons 
in photosynthetic systems. In our model, the influence of the sink is described self-consistently, using a 
non-Hermitian Hamiltonian approach. We include the interactions between the donor and acceptor, and 
between the acceptor and the sink. The classical noise acts only between the donor and acceptor. We 
analyze the dynamics and the efficiency of the ET as a function of the energy gap between donor and 
acceptor, the strength of interaction with continuum, and the noise parameters. We calculate explicitly the 
ET rate and efficiency as a function of parameters. We demonstrate the regimes in which noise assists the 
ET efficiency (in particular, in which the influence of noise significantly increases the efficiency of the ET 
from the "donor-acceptor" subsystem to the sink). 

Our paper is organized as follows. In Section II, using the Feshbach projection method, we introduce 
an effective non-Hermitian Hamiltonian to describe the RC consisting of the donor and acceptor coupled 
to the sink. In Section III, we study the dynamics of the electron transfer without noise. In Section IV, 
we study the decoherence effects caused by the classical noise on the ET efficiency. In the Section V, we 
discuss the obtained results. In the Appendices some important formulae are presented. 
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2 Model description 

We consider a model ("building block" of the LHCs-RCs) of the RC with three sites (protein pigments): 
the first site, \d), is the electron donor (with the energy Ed ), the second site, \a), is the electron acceptor 
(with the energy E a ), and the third site is a "sink", with a continuous spectrum. We assume that the 
acceptor is coupled to the sink, which we first model by a large number of discrete and nearly degenerate 
energy levels, N 3> 1 (Fig. [TJ. The Hamiltonian of this system can be written as, 

v v N 

V .7 cm 



\ d ) T\ 



- Sink \m) 



Fig. 1 A reaction center consisting of the donor and acceptor discrete energy levels, with the acceptor coupled to a 
sink reservoir with a continuous spectrum. 

H t = E d \d)(d\ + E a \a)(a\ + j(\d)(a\ + \a)(d\) 

N N 

+ J2En\n)(n\ + J2 {Vam\a)(m\ + V ma \m) (a\) . (1) 

n—l m—1 

The total Hilbert space can be divided into two orthogonal subspaces generated by two projection oper- 
ators, P = \d)(d\ + \a) (a\ and Q = (|n) (n\), where the P -space is associated with the donor-acceptor 
levels and the Q-space is associated with the sink. These projection operators have the following properties: 
P + Q = 1, P 2 = P, Q 2 = Q and PQ = QP = 0. Then, using the Feshbach projection method Il25ij28l . 
we obtain the effective non-Hermitian Hamiltonian that describes only the "donor-acceptor" subsystem, 

H = E d \d)(d\ + (E a + A(E) - l -T a (E))\a)(a\ + ^(\d)(a\ + \a)(d\), (2) 

where 

i(E) _i r< ,( E , = E J|L. (3, 

n 

To proceed further, we assume that the sink is sufficiently dense, so that one can perform an integration 
instead of a summation. Then we have, 



A(E)- l -T a (E) = J 



V an \ 2 g(E n )dE n 



E - E n +iS 

where g(E n ) is the density of states of the sink. One can show that l29l 

\V an \ 2 g(E n )dE n 



A{E)=V J 



(5) 



E — E n 

T a (E) = 2nJ \V an \ 2 g(E n )S(E - E n )dE n , (6) 

where V denotes the principal value of the integral. 

The exact dynamical evolution of the whole quantum system (RC) is described by the Schrodinger 
equation (we set h = 1), 

^ = HMt)- (7) 



Copyright line will be provided by the publisher 



4 



Alexander I. Nesterov, Gennady P. Berman, and Alan R. Bishop: Non-Hermitian approach 



We assume that at t = the system is populated in the P-space. If the Q-space represents a smooth 
continuum (which is assumed below) one can neglect the dependence of A(E) and T a (E) on E. Denoting 
these functions as A and T a , one can find that the dynamics of the donor-acceptor (intrinsic) states can be 
described by the following Schrodinger equation with the effective non-Hermitian Hamiltonian, li : 

i?^ = WM (8) 

where ip p (t) = Pip(t). Further, it is convenient to rewrite T-L as Ti = H — iW, where 

H = e d \d)(d\+£ a \a)(a\ + j(\d)(a\ + \a)(d\) (9) 

is the dressed donor-acceptor Hamiltonian, W = (l/2)r a |a)(a|, with Sd = Ed and e a = E a + A. 

We define p t (i) to be the density matrix that satisfies the conventional equation of motion with the total 
Hamiltonian, H t : ipt = [H t ,p\. Next, we introduce the projected density matrix as p(t) = Ppt{t)P. 
Then, one can show that p(t) satisfies the Liouville equation, 

ip = [H,p}-i{W,p}, (10) 

where {W, p} = Wp + pW. 

Assume now that the quantum system under consideration interacts with the environment. We use the 
reduced density matrix approach to describe this interaction. To include into the description of the system 
both processes of decoherence and tunneling to the continuum, we introduce the following generalized 
master equation, ip = [H,p] + Cp — i{W,p}, where H is the dressed Hamiltonian, and the Lindblad 
operator, C, describes the coupling to the environment. The commutator of the density operator, p, with the 
Hamiltonian, H, is the coherent part of evolution, and the remaining part corresponds to the decoherence 
process causes by the interaction with the environment. 



3 Tunneling to the sink 



We consider here the quantum dynamics of the ET from the donor |g?) (|1)) to the acceptor \a) (|2)) coupled 
to the sink. We assume that the acceptor is coupled to the A^-level sink reservoir and that the corresponding 
Hilbert subspace is dense and smooth. For description of the tunneling from the acceptor to the sink we 
use the Feshbach projection method described above. This yields the following effective non-Hermitian 
Hamiltonian: 



n = 





' 1 










n 





e + iT 
V 



V 

-e - iT 



(11) 



where Ao = e i + £2 — iT, e — E\ — £2 (s n is the renormalized energy), T = T a /2, with T a being the 
relaxation rate from the acceptor to the sink. 

Region of parameters. The model involves various parameters, which are only partially known. For 
concreteness of the numerical simulations, our choice of the parameters is based on the data taken for the 
ET through the active pathway in the quinone-type of the photosystem II RC l30l (in the units h = 1): 
e = 60ps _1 and 10ps _1 < V < 40ps~ 1 . The parameter T is varied in the interval: lps -1 < V < 5ps _1 . 
But also other values of parameters, e and V, are used in our numerical simulations. (Note, that the values 
of parameters in energy units can be obtained by multiplying our values by h w 6.58 x 10 _13 meVs. For 
example, e = 60 ps _1 « 40meV.) 

In what follows we assume that initially the quantum system occupies the upper level (donor), pn (0) = 
1 (^22(0) = 0). Then, for the diagonal component of the density matrix the solution of the Liouville 
equation ( [TUt is given by (for details see Appendix A), 



,-rt 



(cos T 



Of 



I cos f sm ■ 



-n 



Of 



(12) 
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where fl = \JV 2 + (e + iT) 2 is the complex Rabi frequency, cos 9 = [e + iT)/Cl, and sin 9 = V/Cl. 




t(ps) t(ps) 



Fig. 2 Left panel: The time dependence of the population of the donor site. Right panel: ET efficiency. The 
parameters are the following: blue line (T = 1 ps _1 , V = lOps -1 ), black line (r = lps -1 , V — 20ps~ 1 ), green 
line (T = 5ps _1 , V = lOps" 1 ), red line (T = 5ps~\ V = 20ps _1 ). In all cases e = 60ps _1 . 

The ET efficiency can be defined as the integrated probability of trapping the electron in the sink 0191311 , 



r,(t) = 2T / p 22 (r)dT. 
Jo 

Setting O = Oi + and performing the integration, we obtain for the ET efficiency, 

e -rt 

v(t) = i- T ^ 2 + n2 ^ {(T 2 + n 2 )(T coshn 2 t + n 2 smhn 2 t) 

-(T 2 - fi 2 )(rcos«it- Oisin^it)). 
This yields the following large-time asymptotic behavior: 

^ 1 (T-n 2 )(T 2 + np MMt 



(13) 



(14) 



(15) 



2r(n 2 + n 2 ) 

The numerical results are presented in Fig. As one can see, for these values of parameters, and without 
the action of noise, the ET efficiency approaches a value close to 1 for relatively large times, t > 150 ps. 

Let us consider now the fiat redox potential, e = 0. From the relation £litt 2 = £l\ it follows that for 
e = there are two possibilities: (i) Vt x = 0, fl 2 = \/r 2 - V 2 (V < T); and (ii) VL 2 = 0, fix = y/V 2 - T 2 
(V > r). Using these results we obtain, 
,-rt 



n 2 



■((r 2 (l - cos Hit) -TsinOi*)), V > T 



-n 



v = r 



-rt 

((r 2 (coshfi 2 t - 1) + ^2(^2 + rsinhOit)), V < T 



(16) 



SI 2 



This yields the following asymptotic behavior for the ET efficiency, r/(t) (Tt 3> 1): 

r 2 



1 2nf 

Tt 



-rt 



,-rt 



v > r 
v = r 



i - J_ e -( r -^)* v" < r 
2n 2 



(17) 



Comparing Eqs. (TT3T > - ( fT7l >. we conclude that the highest ET efficiency is obtained for the flat redox 
potential (e = 0), and V > T. 
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3.1 Quantum evolution in the vicinity of the exceptional point 

For the Hermitian Hamiltonian, the coalescence of eigenvalues results in different eigenvectors and the 
related degeneracy, referred to as a "conical intersection", is known also as a "diabolic point" l32l . How- 
ever, in a quantum mechanical system governed by a non-Hermitian Hamiltonian merging not only of 
eigenvalues of the Hamiltonian but also of the associated eigenvectors can occur. In this case, the point 
of coalescence is called an "exceptional point" (EP). At the EP, the eigenvectors merge, forming a Jordan 
block. (For a review and references, see, e.g., 11331 .) 

In the effective two-level system under consideration, the EP is defined by equation £1 = 0. This yields 
e = and V 2 — T 2 = 0. To study tunneling to the sink near a degeneracy, we assume the flat dressed 
redox potential, e = 0. Then, there are two different regimes of the ET depending on the relative values of 
V and r. For V > T, we have a coherent tunneling process (with oscillating probabilities, see Fig. O, 




tips) t(ps) 

Fig. 3 Time dependence of site populations in the vicinity of the EP (red line) for the flat redox potential (e = 0). The 
parameters are chosen as the following: blue line (T = 1 ps _1 , V = 5ps _1 ), black line (r = lps -1 , V = 10ps _1 ), 
green line (r = 5ps _1 , V = 10ps -1 ), red line (r = 5ps _1 , V = Sps" 1 ) corresponds to the exceptional point. 




t(ps) 



Fig. 4 The ET efficiency in the vicinity of the EP (r = Sps" 1 ). Black line (V = 20ps~\ e = Ops" 1 ), blue line 
(V = 10ps"\ £ = 0). Red line: tunneling at the EP ( V = 5ps~\ e = 0). Black dashed line ( V = 2.5ps~\ 
£ = 20ps~ ). Green dashed line ( V = 2.5ps , e = 0). Orange dashed line (V = 2.5ps~ , e = 10ps -1 ). 



n t r . iv\2 ^ Tt v 2 . 2 n t 

pn=e ( cos— + Q^ sm —J ' P22 = e sin — , 



where Sl = \V 2 — T 2 ] 1 / 2 denotes the Rabi frequency. 

On the other hand, for V < T, the tunneling becomes incoherent, without probability oscillations, 

_ vt ( n Q t r n t\2 _ rt v 2 . ^ 2 n Q t 

pn = e r, (rah^|-siiy), p 22 = e rt ^ smh -y- ■ ( 19 > 
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At the EP, Slo = 0, and both regimes coincide. (See Fig. [3] red curve.) In this case, we have the 
following solutions for the probabilities, 

M*) = e- r *(l + ^) 2 , M*) = e" rt (y) 2 . (20) 

The results of numerical simulations of the ET efficiency in the vicinity of the EP are shown in Fig. [4] 
One can see that, for the chosen parameters, the ET efficiency can approach a value close to 1 for short 
times, ~ 2 ps. Note that the coherent tunneling regime (V > T) is more effective for approaching a high 
ET efficiency for short times. (See Fig. [4] black and blue curves.) 



4 Noise-assisted electron transfer 

In this section, we consider ET from the donor, 1 1), to the acceptor, 1 2) , coupled to the sink, in the presence 
of classical noise. Then, the effective non-Hermitian Hamiltonian ( fTTT i takes the form 



U = ^s n \n){n\+^\ rnn {t))\m)(n\ + - ^ \m){n\ - *T|2>(2|, m,n 



1,2, 



(21) 



where A TO „(t)) describes the noise. In our approach, we use a spin-fiuctuator model of noise with the 
number of fluctuators, Af ^> 1 II14II34II351 , The diagonal matrix elements of noise, A„„, are responsible 
for decoherence, and the off-diagonal matrix elements, \ mn (m ^ n), lead to the relaxation processes. 

The approximate equations of motion for the average diagonal components of the density matrix are 
given by (for details see Appendix B) 



J t { Pll {t)) = -W(t)«Pll(t)) - (P22{t)))+0{\V\% 

j t {p 22 (t)) = m((pn(t)) (p2 2 (*)» - 2T{p 22 (t)) + 0(\V\% 
where the average ( ) is taken over the random process describing noise, and 



«K(t) 



-r(t-t'). 



{V{t)V{t')) + {V{t')V{t)))dt'. 



(22) 
(23) 

(24) 



The model of noise. In the following, we restrict ourselves to consider only diagonal noise effects, 
assuming that the noisy environment is the same for the donor and acceptor sites (collective noise). Then, 
one can write \\{t) = gi£(t) and X 2 (t) = g 2 £(t), where £(t) is a random variable describing the stationary 
noise with the correlation function, %(i — tf) = (t;(t)£(t')), and gi ;2 are the interaction constants. We 
describe the noise by a spin-fiuctuator model with the number of fluctuators, Af 3> 1, with the correlation 
function, x( T )' given by lfT4l 



X (r) = a 2 A(E 1 {2 lm T) - £i(2 7c t) 

where E n (z) denotes the Exponential integral 
density of the noise, defined as 

oo 

S(uj) = — x( T ) cos(wr)dr, 

7T J 

o 

is given in |[T4l . 



T = \t-t'\, 

, A = l/ln(7 c /7 m ) and x (0) 



(25) 

a 2 . The spectral 



(26) 



S(u) 



7rwln(7 c /7 m ) 



( arctan ^ 



2 7 „ 



arctan 

V2 7( 



))■ 



(27) 
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This yields the following asymptotic behavior of S(ui): 

2 

a- 



S(u) 



27T7 m ln(7 c /7 

■in J 



1 7m 

7c 



2wln(7 c /7 m )' 

2cr 2 7 c (l - 7 m /7e) 
ituj 2 ln(7 c /7 m ) 



cj <C 27 m , 



2 7m < lj < 2 7c 



w > 27 c 



(28) 



where 7,„ and j c (7 m <C 7 C ) indicate the boundaries of the switching rates in the ensemble of random 
fiuctuators. As one can see from Eqs. (f28}, for u> <C 27 m the spectral density of noise, S(uj), describes the 
white noise. In the interval of frequencies, 2j m <C u> <C 27 c , one has the 1// noise: 5 ~ 1// (/ = us/2ir). 
And for w ^> 27 c , we obtain the Lorentzian spectrum. 

Choice of parameters. The correlation function includes, besides the amplitude, a, two fitting parame- 
ters: 7 m and j c . Taking into account available theoretical and experimental data (37—39), we have chosen 



10- 4 ps _1 , 2 7c 



lps . Note that as our 



Is practically does not change 



in our numerical simulations the following parameters: 27^ 
results demonstrate, a decrease of the left boundary, j m , even up to 7^ 
ET rates. We also introduce the notation: D = \gi — (72 1 ■ 

The spectral density of noise corresponding to Eq. (l27t and its asymptotic behavior given by Eq. 
is presented in Fig. [5] 



S{f) 




In S(f) 




0.10 0.15 
THz 



0.00001 0.00003 0.00005 0.00007 
THz 



Fig. 5 The spectral density of noise a — 1. Left panel: S(f) given by exact formula, Eq. &27\ (blue line), and the 
asymptotic formula d28t : S(f) ~ 1 // (red line). Right panel: In S(f) as a function of frequency. 



4. 1 Influence of noise on the electron transfer rates (T = 0) 
For r = 0, we obtain the exact solution of the system (1221 and (l23l . 

(pn(t)> = 5 + (pn(O)-0e- 2 /o t W j (29) 

(m{t)) - \ + (p 22 (0) - e - 2 /o«(*')*'. (30) 

It follows that, independent of the initial conditions and the nature of noise (producing decoherence or 
relaxation), in the limit, t —> oo, the presence of noise produces equal populations in the two levels. (See 
Fig. |6l left panel.) 

The computation of the ET rate, 9t(i), yields, 

£Ft(t) = ^ / e l£T <e 8K(r) )dr, k(t) = -D / (31) 
4 j-t Jo 
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To proceed further, we use the first order cumulant expansion (the Gaussian approximation) to evaluate the 
generating functional (e lK ^Y The computation gives 

<e IK ) = e -<- 2 >/ 2 = cxp ( - D 2 J dt' j dt" X (t' - t")) ■ (32) 

Let us assume that initially the system occupies only the upper level (donor), pn(0) = 1. Then, if 
,9i = .92 the solution for the diagonal components of the density matrix takes the form, 

/ / ^ 1 1 I n V 2 . 2 et\ , , 1 1 / V 2 . 2 et\ 
<PnW) = 2 + 2 CXp (" V 2")' (-° 22 W) = 2 _ 2 CXP (" 2 72" S 2")- (33) 

One can see that up to the first order in the dimensionless parameter, V 2 /e 2 , the approximate solution ([33} 
coincides with the exact solution given by Eq. (fT2l (with V = 0). In this case, the effect of the collective 
noise vanishes. 




0.5 1 1.5 2 0.05 0.10 0.15 0.20 0.25 

tips) t(ps) 



Fig. 6 Left panel: Effects of the diagonal noise on the time dependence of the site populations: pn(t) (blue line) 
and p22(t) (red line). Right panel: Blue line describes the time dependence of the ET rate, SR(f), given by Eq. ( 13 1 1 . 
2jm = 10 _4 ps _1 . Green dashed line corresponds to 2j m = Is -1 . Red line corresponds to the asymptotic formula 
d35t . The parameters are: Da — e — 60ps _1 , 2j c = lps -1 and V = 20ps~ 1 . 



Relation to the Marcus' theory. The asymptotic ET rate for L = is defined as 1Z = lim t _ i . 00 *R(t). 
Using Eq. (IBTb . we obtain, 



V 2 f°° 

K=— dtexp(ist-Q(t)), 

^ J — OO 



(34) 



where Q(t) = D 2 J* dt' J* dt"x(t' - t"). To evaluate 6(i), we use the approximation, x(t) « x(0). 
(Note that x(0) = a 2 .) This yields Q(t) )c (aDt) 2 jl. Performing the integration over t in (l34l , we obtain 



V 2 HhT / e 2 



K = TV5v cxp ^wj- (35) 

In Fig. [6] we compare the results of numerical calculations (blue line) of relaxation rate, £R(t) described 
by by Eq. d3TT > with the asymptotic formula (f35T >. One can see a good agreement of the asymptotic rate 
defined by Eq. d35l ) with the formula (|3"T1 i. 

In the case in which the number of thermally excited fluctuators, Mt 3> L the dispersion, c 2 , is a linear 
function of the temperature, so that a 2 = P$kT [35). Inserting this expression into Eq.(l35l>. we obtain 
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where A = D 2 Pq/2 and \V\2\ — V/2, Comparing our result with the Marcus formula [|8], one can see that 
the classical noise results in a large-time asymptotic of the ET rate which can be expressed in the form of 
the Marcus-type formula. 

4.2 Noise-assisted electron transfer in the reaction center (T 7^ 0) 

We consider here noise-assisted ET to the sink described by the equations of motion (1221 and (1231 for the 
averaged components of the density matrix: 



j t (p n (t)) = -X(t)(( Pll (t)) - (p 22 (t))), 
j t {p 22 (t)) = m((pn(t)) (p 22 (t))) 2T{p 22 (t)), 



(37) 
(38) 



where SK(f) = (1/2) V 2 /„ e" r * coser (e iK ^)dr. Further, we use the Gaussian approximation to evaluate 
the generating functional: 



(e-«) = exp ( - ( 9l - g 2 ) 2 ^ dt> £ dt" X (t> t") 



(39) 



Performing the integration over t in Eq. t24l . we obtain for T 7^ a generalization of asymptotic expression 
for 1Z given by Eq.(l35ll for finite T, TZr = limt_j.oo 9^(f ), as 

where erfc(z) denotes the complementary error function, crfc(z) = 1 — erf (z) 13 61 . The dependence of 
1Zy as function of the amplitude of noise, Da, is presented in Fig. [7] 




20 40 60 80 100 120 140 160 180 200 
Da 



Fig. 7 The function TZr vs. the amplitude of noise, Do (e = 60ps~\ V = 20ps _1 ). Black line (T = 10ps _1 ), 
green line (T = 5ps _1 ), blue line (T = lps -1 ), red line (r = Ops -1 ). 



Using these results and taking the initial conditions as pii(0) = 1, we obtain the approximate solution 
ofEqs. ||37) and d38]) 



-75.it 



2^K 2 +T 2 , 

(P 22 (t)) 



2 2^U 2 +T 2 



-n 2 t 



e -n 2 t _ e -Kit 



where 72-1,2 = 7£r + T ± + T 2 . Inserting (l42l into Eq. (fT~3b . to we obtain for the ET efficiency 

2 



17(f) = 1 - e 



[ cosh + sinh 

2 TZi-TZ 2 2 



(41) 
(42) 

(43) 
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This yields the following asymptotic behavior of the ET efficiency 

As can be seen from Eq. (l43l > there are two ET rates, IZi and 1Z2- The asymptotic behavior of ET efficiency, 
r](t), is defined by the lowest ET rate, 7^2- 




t{ps) t(ps) 



Fig. 8 The time dependence of the site population (left panel) and the ET efficiency (right panel) in the presence of 
noise. Left panel: (i) V = 20ps -1 , blue and black lines correspond to, pn(i) and p22{t), respectively; (ii) V = 
40ps _1 , red line corresponds to pn(t), and green line corresponds to p22(i). Right panel: The ET efficiency, r/(t), in 
the presence of noise. Green line: (V = 20ps _1 , Da — 30ps _1 ); red line: (V = 20ps~ 1 , Da « 60ps _1 ); black line: 
(V = 20ps~\ Dcr = 120ps _1 ); blue line: (V = lOps" \ Da = 60ps _1 ). In all cases: e = 60ps"\r = lps -1 , 
2j m — 10 _4 ps _1 and 2j c = lps -1 . 

In Fig. [8] we present the results of numerical simulations for V = 1 ps _1 , and for a given sharp redox 
potential, e = 60 ps _1 . Red line corresponds to the amplitude of noise which was optimized by using the 
modified Marcus-type formula given by Eq. (l40l . (See Fig. [7]) As one can see from a comparison of the 
results presented in Fig. |2](for a sharp redox potential) and Fig. [HJ influence of noise with amplitudes near 
to optimal value significantly accelerates the ET to the sink. 

5 Conclusions 

In this paper, we model quantum electron transfer dynamics in a photosynthetic reaction center consisting 
of three elementary pigment units. Two of them, a donor and an acceptor, are represented by localized 
sites of protein pigments with discrete energy levels. The donor interacts with the acceptor through the 
corresponding matrix element. The third protein pigment (sink) has a continuous energy spectrum, and is 
described by two parameters: its density of states and its strength of interaction with the acceptor. The sink 
is described self-consistently, by using the Feshbach projection method on the "donor-acceptor" intrinsic 
states, within a non-Hermitian Hamiltonian approach. We apply our results to the quantum dynamics of 
the electron transfer in the active branch of the quinone-type PSII reaction center. 

The collective external noise produced by the environment of the proteins acts on the "donor-acceptor" 
sub-system. Usually, the presence of noise acts as an incoherent pump in the system under consideration. 
But, as our results demonstrate, the simultaneous influence of both noise and the sink, significantly assist, 
under appropriate conditions, the quantum efficiency of the electron transfer. We derived the expression 
for electron transfer rate which describes the tunneling to the sink, in the presence of noise. We calculate 
explicitly the corresponding region of parameters of the noise assisted quantum electron transfer for sharp 
and flat redox potentials, and for noise described by an ensemble of two-level fluctuators. 

Our results show that even in this simplified model, the quantum dynamics of the electron transfer to 
the sink can be rather complicated, and depends on many parameters. Further analytical research and 
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numerical simulations are required to extend our approach for (i) complicated dependencies of the density 
of states on energy in the sinks for fiat and sharp redox potentials, in the presence of noise and thermal 
environments, and (ii) more complicated LHCs-RCs complexes. The problem of the electron transfer 
optimization also requires further analysis. 
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A Solution of the Liouville equation 

The solution of the the Liouville equation ([Toi l, ip = [H, p] — i{W, p], is given by M401I411 . 

flt\ fit 
. — (7i — i sin 9 sin — G2 
2 J 2 

nt\ n . fit 

- 1 G2 — 1 sin 9 sin — Gi 



Pu{t) 

P22(t) 



,-rt 



/ fit 

(cos T 



fit 
Y 

fit 



l cos f Sill 



2 / 



(45) 
(46) 



-17 



fit 



nt 



1 cos t) sm ■ 



1 cos f sm 



P2lW= e - rt (( 



2 

fit 
~2 I 
fit 



Gi — z sin & sin— G2 



• • a ■ nt n 
2 — J sm a sin — C\ 



fit , 

cos 1 cos t> sm — 

2 2 



C'i — i sin 8 sin 



fit 



fit 



2 

fit 

T 
fit 



1 cos V sm 



2 7 



~] C2 — i sin 9 sin 



(47) 



(48) 



where cos 9 = (e + iT) /fl, sin 9 = V/fl, and fl = \JV 2 + (e + iT) 2 is the complex Rabi frequency. The 
constants, C\ and C2, are defined from the initial conditions as follows: pn(0) = |Gi| 2 , ^22(0) = | C2 1 2 , 
Pi 2 (0) = C*C 2 and p 2 i(0) = C*d. 

Let us assume that initially only the acceptor site is occupied, so that pn{0) = and ,022(0) = 1. This 
yields, 



-17 



fit 



sm v sm ■ 



Pi2{t) =ie~ rt 
P2i(t) = - 



rt 



P22{t) 



sm a sm 



sm a sm 



-VI 



t fit 
(cos T 



fit\* 


/ ot 




[ C0S Y 


flt\ 


( fit 




l cos T- 


T • 





+ i cos sm — 



1 cos f sm 



Of 
> 

flt\ 
2j' 



fit 



i cos sm ■ 



(49) 
(50) 
(51) 



Taking V = 0, we find pn(t) = Pi 2 (t) = p2i(t) = 0, and for P22(t) we obtain the Weisskopf-Wigner 
formula for an irreversible decay, P22 (t) 



-2Tt 



Presenting fl = fli + ifl2 = \/p + iq, where p = V +e —T* and q = 2eT, we obtain Oj — fl% 



and flifl2 = q. This yields, fix = ±^^Jp+ ^Jp 2 + q 2 , O2 = — p+ \Jp 2 + q 2 , where the 



V2 
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upper sign corresponds to p > 0, and the lower sign corresponds to p < 0. Using these results, we obtain 
for /022(f) the simple analytical expression, 

V 2 e~ rt 

/022(f) = 2 ^2 + ^2) (cosh» 2 f - cosfiif). (52) 

B Equation of motion for the average density matrix 

In the interaction representation, considering the off-diagonal elements as perturbations, so that % = 

U a + V{t) - iW, where 

H = J2 £ n\n)(n\+J2 x nn(t)\n}{n\, (53) 

n n 

V(t)=J2( V mn + Xmn(t))\m){n\, W = r|2)<2|, (54) 

we obtain the following equations of motion, 

Pn = Kh^i - V12P21), P22 = i(hiVi2 - V21P12) - 2Tp 2 2, (55) 
p 12 = iVi 2 (pii - P22) - Tp 12 , P21 = iV 2 i(pu - P22) - Tp 2 i, (56) 

where p = T{e i ^o H °^ dT )pT(e~ i ff °( T ) dr ) and V = T(e l ti H ^ T ) dr )VT{e~ i ^o H (T)dry 
Using Eqs. 05] )- ( l56l l. we obtain 

pu(f) = pii(0) + i f (/5i 2 (f')V2i(f') - X> 12 (f')/52i(f'))rff', (57) 
Jo 

/022(f) = MO) + * / e" 2r(t " t ' ) (P2i(f')^i2(f') - V2i(f')Pi2(f')), (58) 
Jo 

/012(f) = pa(0) + i f e- r ( t - t ')y 12 (f')(pii(f') - Pi2{t'))dt' , (59) 
Jo 

/021(f) = 021(0) + * I e-^-^VtiWifiixtf) - P22(t'))dt' . (60) 
Jo 

We assume that initially pi 2 {0) = P2i(0) = 0. Now inserting d57] > - (l60l into Eqs. (f55l > - d56l l. and 
taking into account that pn = p\\ and P22 = P22, we obtain the following system of integro-differential 
equations, 



Pu = - e" 
Jo 

P22= /V r(t 







r(t - 4,) (v 21 {t)V l2 {t') + V2i(t')V 12 (t)) (pn(f) - P2 2 (t'))dt', (61) 
(v 2 i(f)V r i 2 (f') + V2i(f')^i 2 (f)) (pn(f') - P22(t'))dt' - 2Tp 22 (t), (62) 



Pi 2 (f) = - / (l + e- 2r ( 4 -*')) (y 21 {t>)p 12 (t>) - U 12 (f')P2i(f')) Vi 2 (*)dt' - I> 12 (f) 
+ iV r i2(f)(pii(0)-/02 2 (0)), (63) 

P2l(f) = - j (l + e- 2 ^ 4 -*')) (>2l(f )/5l2(f) - Ul 2 (f')P2l(f')) V2l(f)rff - I>2l(f) 

+ ^2i(f)(pii(0)-p 22 (0)). (64) 
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(Pai(*)> = - / (l + e- 2r ^)((v 21 (t')p u (t')-V 12 (t')^ 



For the average components of the density matrix this yields, 

J t (Pn(t)) =- j^e- T ^((v 21 {t)V 12 (i!) + V 21 {t')V 12 {t)) (p u (f ) - p 22 (t')))dt' , (65) 

J t (p2 2 (t)) = ^e- r ^((v 21 (t)V w (t') + V 21 (t')V 12 (t)) ( Pll (f) - p 22 {t'j))dt> 
-2T(p 22 (t)), (66) 

|W 2 (*)) = (i + e- 2r ^- t, ))((v r 2 i(f)Pi 2 (*') - fcaWftiCO)^*))^ - iWt)> 

+ «(V r i 2 W)(pii(0)-p 22 (0)), (67) 
d 

+ *<V r 2 1 (t)>(pii(0)-^ 22 (0)) ) (68) 

where the average ( ) is taken over the random process describing noise. Generalization of the obtained 
results for the case |5i 2 (0) 7^ and p 2 i(0) ^ is straightforward. 

In the spin-fluctuator model of noise with the number of fluctuators M ^> 1 one has the following 
relations for the splitting of correlations lfj"4l . 

( (V21 {t)V 12 (O + V 21 {t')V 12 (i)) (p n {t') - p 22 (t'))) = 

(<F 2 i(t)Fi 2 (tO> + <^i(t')Vl 2 (i)>)(<Pu(t')>-<P2 2 (0>) ) (69) 

and so on. Next, using the second-order cumulant expansion, we obtain the following system of differential 
equations for the average components of the density matrix, 

^<Pu(*)> =-f*e- T ( t - t '\(V 21 (t)V 12 (t')) + (V 21 (t')V 12 (t)))dt'(( Pll (t)) (p 22 (t))) 
+ 0(V 4 ), (70) 

dt 

-2T{p 22 {t))+0{V A ), (71) 
d 



-{p 22 (t))= re- r( *- t ' ) ((^iW^i2(i')) + <^i(i')^i2(t)))^«PiiW)-<P2 2 W)) 



-(p 12 (t)) = i(V 12 (t)}( Pll (0) - P22 (0)) ~ (l + e- 2r ^) {V 12 {t)V 21 {t'))dt'{p 12 {t)) 
+ J (l + e- 2T ^)(V 12 {t)V 12 {t'))dt'(p 21 {t)) - T(p 12 (t)) + 0(V 4 ), (72) 
j t (fci(t))=i(V 21 (t))(p u (0)-p 22 (0))- /' t (l + e- !ir (*- 1 '))<V a i(t)y a i(f))d«'<Pi2(t)> 



f ^ + e - 2r (*-* , ))(V ai (t)V 12 (*'))tft'(AJi(«)> - iWt)) + 0(^ 4 )- (73) 
We rewrite the equation of motion for the diagonal components of the density matrix as, 

j t ( Pll (t)) = - (74) 

j t (p 2 S)) = <R(t)((p n (t)) - <P2 2 (0» - 2r(/5 22 (t)), (75) 
where 3t(t) = /* e -r(*-f) «F 21 (t)F 12 (*')) + (V 21 {t')V 12 (t)))dt>. 
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